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The binding energies of baryonic systems with baryon number i? = 2,3 and 4 
possessing heavy flavor, charm or bottom, are estimated within the rigid oscillator 
version of the bound state approach to chiral soliton models. Two tendencies are 
noted: the binding energy increases with increasing mass of the flavor and with 
increasing B. Therefore, the charmed or bottomed baryonic systems have more 
chances to be bound than strange baryonic systems discussed previously. The flavor 
symmetry breaking in decay constants F is considered which is especially important 
Ch , for baryonic systems with bottom quantum number. 
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CN . 1 Introduction 

^ ■ Many efforts have been done lately to investigate the properties of baryonic systems 
\^ I {BS) with nonzero strangeness, first of all the possibility of the existence of states 
0> ■ stable relative to strong decays. 

Recently some of the predictions of theory began to find experimental con- 
I firmation. The near-threshold enhancement in AA system observed in Q] can be 
00 ' interpreted as a component of 27-plet obtained from the bound SU{2) torus-like con- 
^ ! figuration with i? = 2 by means of collective coordinates method described in 0, ^ . 
' Similar enhancement in AA^ system has been observed many years ago in the kaon 
production reaction on nucleons and confirmed also in Ap scattering 0. It can 
Qh! belong to 27-plet or to antidecuplet of dibaryons. The singlet A^A^ scattering state 
with isospin T = 1 belongs to the 27-plet (for review of theoretical predictions in 
B = 2 sector see, e.g. 0). Analogous results are obtained in more conventional 
potential approach as well. 
^ ■ The question if the BS with flavor different from u and d can exist, is more 

general, of course. Charm, bottom or top quantum numbers are also of interest. 
Their consideration can be performed in the framework of chiral soliton models, in 
particular, the bound state approach to heavy flavors proposed in and developed 
in 0-[|lO[- Although charmed and bottomed BS have less chances to play some 
important role in astrophysics than the strange ones (it is not excluded, however!) 
their studies can be very useful for understanding of the peculiarities of nuclear 
matter fragments with unusual properties. It might be similar to heavy quarkonia 
which studies were very important for development and checking of QCD itself. 

Here the baryonic systems with heavy flavors are considered within the rigid 
oscillator version of the bound state approach to strange baryons proposed by Ka- 
plan and Klebanov and used later in This model has deflnite advantages 
before collective coordinates quantization method when heavy flavors are included 

^The work supported in part by the Russian Fund for Fundamental Research, grant 95-02-03868a and 
presented at the Workshop on Science at Japan Hadron Facility (JHF98), KEK, Tsukuba, March 4-7, 1998 



into consideration, first of all, because of its simplicity. However, some apparent 
drawbacks are present also. 



2 Rigid oscillator model 

The ansatz for the chiral fields used in |10[ is: 

U{r,t) = R{t)Uo{r)RHt), R{t) = A{t)S{t), (1) 

where Uo is SU{2) soliton embedded into SU{3) in usual way (into left upper corner), 
A{t) e SU{2) describes SU{2) rotations, S{t) e SU{3) describes rotations in the "charm" 
or "bottom" direction. For definiteness we shall consider the extension of the {u,d) 
SU{2) Skyrme model in charmed direction, when D is the field of £)-mesons. But it 
is clear that quite similar the extension can be made in bottom and top direction. 

S{t)^exp{iD{t)), D{t)^ Da{t)Xa, (2) 

a=4,...7 

Aq are Gell-Mann matrices of iu,d,c) or {u,d,b) SU{3) groups. The {u,d,b) SU{3) 
subgroup is quite analogous to the (u, d, s) one, for (u, d, c) subgroup simple redefiniton 
of hypercharge should be made. D4 = {D° + D")/^/2, D5 = i{D" - D°)/V2, etc. 

After some calculation the well known Lagrangian of the Skyrme model in 
the lowest order in field D takes the form [0, [T0[| : 



L = -M,i^B + AQf^bD^D - TB{m]j - ml)D^D + i^{D^D - D^D). (3) 

Here D is a doublet formed by D° and mesons, and we maintained our former 
notation for the moment of inertia for the rotation into "strange", "charm" or 
"bottom" direction 9c = 6& = 9^ = 9f. This moment of inertia has simple analytical 
form for arbitrary starting SU{2) skyrmion, regardless its symmetry properties: 

Qf,b = ^ /(I - ^f)[P' + ^(.(dfr + s){daf + s}sl{d[(3Y)]d\ (4a) 

Ft, and e are the parameters of the model. The general parametrization of the 
SU{2) skyrmions has been used here, U ^ Cf + SfTfi with = Ca, = SaC/s, Uy = SaSp, 
Sf ^ sinf, Cf = cos/, etc. For the axially symmetrical ansatz (5 = n0, is the azimuthal 
angle, and 6f,b takes the form drawn in [ITT 



e^.B-J y'(l-c/)[F2 + l((/,/) + 4(a,a) + ^44)]rdrdz-^^ + 9|fB, (46) 

(/:/) = {df /drY + {df /dzY, r and z being cylindrical coordinates. The quantity Tb 
defines the contribution of the mass term in the Lagrangian: 

rs = ^ /(i-c/)dV (5) 



Numerical values of 9f,b, Tb and some other quantities are shown in Table 1 below. 



The term in (3) proportional to NcB appears from the Wess-Zumino-Witten 
term in the action and is responsible, within this approach, for the splitting between 
excitation energies of charm and anticharm (flavor and antiflavor in general case) 



[P]-[|10|]- Nc is the number of colors in the underlying QCD, in all other cases here 
the index c means the charm quantum number. B is the baryon number of the 
configuration which can be written in terms of the functions /, a and /3 as 

B = J s)s^{dfdadl3)d\, (6) 

In other words, it is the Wronskian of the system described by 3 profiles, /, a and 
(3 0. For the axially symmetrical configuration possessing also symmetry 
B = n(/(0) - /(cx)))/7r = n for configurations of lowest energy. 



z 



As it was noted in deviations of the field D from the vacuum decrease 
with increasing mass m^i, as well as with increasing number of colors Nc, and the 
method works at any - for bottom and even top quantum number also. The 
zero modes quantum corrections due to rotation with the matrix A{t) have the order 
of magnitude N~'^ and are not crucial but also important (see also section 4). 



3 Flavor excitation frequences 

After the canonical quantization procedure the Hamiltonian of the system takes the 
form: 

Hb = M^^B + TT^^ntn + {Tsfnl + ^§^)D^D - ^^^iD^U - U^D). (7) 



ml The momentum n is canonically conjugate to variable D. Eq. (7) 
describes the oscillator-type motion of the field D in the background formed by 
the (u,d) SU{2) soliton. After the diagonalization which can be done explicitely 
according to |10[ the Hamiltonian can be written as 

Hb = McLB + 0JF,Ba^a + uip.Bb'^h + 0{l/Nc) (8) 

with a\ 6^ being the operators of creation of charm and anticharm (bottom and 
antibottom) quantum number, ujfm and u)f,b being the frequences of heavy fiavor 
(antifiavor) excitation. D and n are connected with a and h in the following way 
1,0: 

D'= {a' + h% n'^ ^^f^"'" (a'-6t') (9) 

y'NcBnF,B 2i 

with 

f^F.B = (1 + mmlTBeF,B/iN,Bf)'/^. 
For the lowest states the values of D are small, 

D ~ [iTBBF.Bml + NlB'^I^Y^'\ 
The flavor (antiflavor) excitation frequences w and uj are: 

NcB , . _ NcB , , , , 

(^F,B = -^ (mf,b-1), '^■F-B = -J?; + (10) 

SUfs o'C>f,B 



It should be noted that the difference lufm ~ ^f,b = NcB/{AQfm) coincides in the 
leading order in Nc with that obtained in the collective coordinates approach [jl2|, [13 



Indeed, in the collective coordinates approach the zero-modes energy of the soliton 
rotated in the SU{i) configuration space and depending on the "fiavor" inertia 9^,^ 
can be written as: 



Erot{'S>F,B) 



46 



F,B 



■ [N,B + nqq{N,B + 2nqq + 2 - 2Tr)] 



(11) 



where n,^ is the number of additional quark-antiquark pairs present in the quantized 
state, NcB + inqq = p + 2q, p,q are the numbers of indices in the spinor describing the 



qq)/2 is the so called right isospin characterizing irrep (see [p? 

and 



121 where (11) was obtained for 



qNcB in (11) coincides with the difference of 



SU{3) irrep, Tr = {p- 

where the B = l,nqq = case was considered 
Nc = 3). The term proportional to 
LUFM " ^F,B in (10). 

For the difference of the frequences of excitation in cases of S > 2 and 5 = 1 
systems we obtain: 



Auj 



tuf 



6 



1/2 



6 



F,B 



1/2-1 



(12) 



It is proportional to the heavy quark mass and is positive if ri/6_F,i > Ts/e^.s- 
For B = 2, 3 it is really so. The characteristics of S'U{2) toroidal solitons with baryon 
numbers B = 2,3 and 4 have been calculated previously [jl^. For B = 2 they coincide 
with good accuracy with those given later in [|10[ . For greater baryon numbers some 
configurations of lower energy have been found [|T5| , p!6| , but necessary quantities 
like eF,s and are absent, still. 

As a result, the binding energy of heavy fiavored dibaryons, tribaryons, etc. 
increases in comparison with strange fiavor case, as it can be seen from the results 
of numerical estimates shown in Table 1. 



4 zero modes corrections and binding energies estimates 



The 1/A^c zero modes quantum correction to the energies of BS can be estimated 
according to the expression [^, p!0[| : 



26 



T.B 



■ [cF,BTr{Tr + 1) + (1 - cf,b)I{I + 1) + {cf,b - cf,b)T{T + 1)] , (13) 



where / is the isospin of the BS, Tr is the quantity analogous to the "right" isospin 
Tr in the collective coordinates approach 0, |ll|. 



CPS = 1 



6 



t,b 
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F.BlJ'F.B 



-{fJ-F,B - 1), 



CF,B = 1 - 



and fr = + f . 
6t,s 



Qf,b{p^f,b)' 



{fiF,B - 1) 



(14) 



In the rigid oscillator model the states predicted are not identified with definite 
SU{3) or SU{A) representations. However, it can be done, as it was shown in []T^ . 
The quantization condition {p + 2q)/3 ==5 for arbitrary Nc is changed to {p + 2q) = 
NcB + 3nqq. For example, the state with c = 2, 1 = and = should belong to the 
27- plet of {u, d, c) SU{3) group, if Nc = 3, see also [|0l- For 27-plet of dibaryons Tr = 1, 



for antidecuplet Tr = 0. For 35-plet of tribaryons Tr = 1/2, for arbitrary {p,q) irrep 
which the BS befongs to Tr = p/2 if n,, = 0. / and T take the fowest possible values, 
or 1/2 in our case. If 9^ oo Eq. (13) goes over into the expression obtained for 
axially symmetrical BS in collective coordinate approach []ll| , in realistic case with 
Qt/^f — 2.7 the structure of (13) is more complicated. 

The quantum correction due to usual space rotations, also of the order of 1/Nc 
is exactly of the same form as obtained in ||11|| , 



see 



The binding energies 
shown in Table 1 are defined relative to the decay into B baryons, nucleons or 
flavored hyperons. The binding energy, e.g. of B = 4 state relative to 2 dibaryons 
will be smaller or negative. Since we are interested in the lowest energy states we 
discuss here the baryonic systems with the lowest allowed angular momentum, J = 
for B = 2, 4, and J = 3/2 for B = 3. The latter value is due to the constraint because 
of symmetry properties of the configuration. The value J = 1/2 is allowed for the 
configuration found in ||15|| . 

For i? = 3 and 4 toroidal configurations we used here do not correspond to the 
minimum of static energy, but only for such configurations the necessary quantities, 
Qf,b, are known. For 5 = 3 the toroidal configuration does not differ much in 
energy from the tetrahedral one which is known to be the configuration of minimal 
energy 



(The masses of stranglets obtained from bound skyrmions with B 
up to 17 [rofhave been estimated recently in [|T^ in the bound state soliton model.) 
For 5 = 4 the difference is large, ~ SOOMev in energy. However, it would be incorrect 
to decrease all B = 4 energies by SOOMev and increase the binding energies, because 
other characteristics of solitons and, therefore, the excitation energies Wc and ujb also 
change. Some reasonable extrapolation for i? = 4 is shown in Table 1. 
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UJc 




es=-2 




ec=2 
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efc=-2 
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0.865 


1.86 


5.14 


5.14 


3.98 


0.200 


1.18 


3.66 












2 


1.656 


3.79 


10.55 


16.45 


7.80 


0.196 


1.15 


3.62 


0.096 


0.16 


0.15 


0.18 


0.19 


3 


2.523 


6.16 


16.85 


37.85 


12.85 


0.205 


1.17 


3.63 


0.12 


0.22 


0.23 


0.27 


0.27 


4 


3.446 


8.84 


23.65 


72.5 


18.80 


0.215 


1.19 


3.68 


0.18 


0.23 


0.21 


0.25 


0.25 


4* 


3.140 










0.196 


1.15 


3.62 


0.52 


0.58 


0.61 


0.60 


0.65 



Table 1. The static characteristics of the B = 1 hedgehog and toroidal solitons with 
B — 2,3,4 Q: Mci^B in Gev, moments of inertia 6f,b = &c = ©6, ©t, ©j and T in Gev~^. 
The excitation frequences Ws,c,6 - in Gev. The binding energies (in Gev) of baryonic systems with 
B = 2,3,4, s = —2 {es=-2), charm c = 1,2 (ec=i,2) and bottom b = —1,-2 {eb=-i.-2) are shown. 
The parameters of the model = 108Afew, e — 4.84 [||. The line B = 4* shows the binding energies 
for B — 4 configuration found in fisl , ^ with extrapolation ljb=4 = ^b=2- The uncertainty of these 
estimates within our choice of the model and configurations is ~ 0.02G'eti. 

The binding energy of the deuteron-like state within this approach is 0.16 
Gev, the binding energy of the NN scattering atate with isospin T = 1 is 0.127 Gev. 
If we assume that the nonzero modes quantum corrections - due to vibration and 
breathing modes, as well as 1-loop corrections - are approximately the same for all 
states with the same B then we should renormalize the energy of each state ob- 
tained in this way, i.e. we should subtract just 0.127G'eu from the binding energies 
shown in Table 1 (the energy of the virtual level, 0.067 Mev, is negligible, of course). 



Renormalization of this type was done also previously in [|TT[. It is clear that after 



this renormalization the state with strangeness = -2 becomes unbound, 0.03 Gev 
above the threshold, but the states with charm or bottom remain bound. These 
estimates are crude, of course, because the binding energy of the deuteron is about 
0.03 Gev within the same approach. However, they show clearly that BS with 
charm and bottom have more chances to be bound relative to strong interactions 
than strange BS. 



5 Symmetry breaking in flavor decay constants 

The estimates of binding energy were made in ||^, |10[ with the constant Fx = Ft,, and 
similar assumption was made above to get the results shown in Table 1. However, 
symmetry breaking takes place not only in the masses of baryons, but also in decay 
constants, since generally Fp is different from F^. As it was shown by Riska and 
Scoccola [|TB[ , the large enough values of Fp allow to remove the apparent overbind- 



ing of heavy mesons by SU{2) solitons which is characteristic for the bound state 
models with Fp = F^,. 

For kaons, Fk = l.22F^, and this is important for the description of mass 
differences inside the octet and decuplet of baryons. Therefore, it seems necessary 
to investigate the binding energies of BS also for the case of this kind of flavor 
symmetry breaking [FSB). 

The mass term in the Lagrangian (3) should be changed to 

LM^-{mlF^/F^-ml)rB. (15) 

It is assumed that the pion mass term is included into the classical mass of the 
soliton Mci- Corresponding change should be made in the Hamiltonian (7) also. The 
second order term in the Lagrangian should be modified also, but the corresponding 
contribution to the mass or energy is much smaller than that given by (15) because it 
is proportional to msf in the integrand and cancellations take place in the integration 
over tfr, see 0, ^ 



It is not difficult to calculate the "fiavor" moment of inertia in the case of 
FSB: 

&F,B ^ eF,B + [Fl/F^ - l)rB/4 (16) 

The modified inertia Qp can be calculated easily according to (16) for any Fp since 
9;- and are known. 

The following expression can be obtained for the excitation energy: 

rnp 

Here rp = Fp/Ft,, Qp^B is the Skyrme term contribution to the moment of inertia, so 
we can write: 

eF.i3 = ef^B + r|rB/4. (is) 

For our choice of the model parameters Off^B is a bit smaller than rs/4. 
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2.346 


3.104 


4.845 


0.255 


1.466 


4.708 


0.480 




0.883 




0.980 




2 


4.742 


6.228 


9.640 


0.250 


1.447 


4.671 


0.472 


0.11 


0.880 


0.16 


0.982 


0.20 


3 


7.729 


10.18 


15.80 


0.260 


1.464 


4.684 


0.498 


0.15 


0.892 


0.22 


0.984 


0.28 


4 


11.13 


14.71 


22.94 


0.270 


1.481 


4.713 


0.528 


0.20 


0.903 


0.26 


0.986 


0.30 


4* 








0.250 


1.447 


4.671 




0.52 
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Table 2. The flavor inertia Qf,b (in Gev) for Fk/F^ = 1.22, Fc/F^ = 1.5, Fb/F^ = 2. 
Excitation frequences and some binding energies for |F| = 2 baryonic systems with flavor symmetry 
breaking in the constants Fp - in Gev. The quantities Cf,b defining the suppression of zero-modes 
quantum corrections are shown also. 

We took here rc = Fd/F^, = 1.5 and n = Fb/F^ = 2 because these values of decay 
constants ratios aUow to get the masses of lowest baryons Ac and not very far 
from the observed values, rc does not contradict to the experimental restriction. As 
it can be seen from Table 2 the flavor symmetry violation in meson decay constants 
leads to moderate increase of the binding energies of BS with different flavors, and 
the tendency of increase of with increasing mass of flavor becomes more striking. 

It is of interest to look at the case of large FSB, when the ratio rp = Ff/Ft, > 1. 
In this case it is easy to flnd out what is the effective potential for the heavy meson 
in the SU{2) soliton background. The following expression can be obtained: 

UJF,B ^mF[l- 2e^fi5/(rBr|,)] - ml/{2mFrl) - N^B/{2TBr%) (19) 

This defines the average potential of the heavy meson bound by SU{2) soliton. Ob- 
viously, iOF,B rnp when rp oo. 



6 Conclusions 

To conclude, we estimated the binding energies of dibaryons, tribaryons and tctra- 
baryons with nonzero charm and bottom. For the top quantum number similar 
results can be obtained, but the spectroscopy of mesons, baryons and baryonic sys- 
tems with t-number will not be available, probably, because of the large width of 
the top quark. 

For rt = Ft /Ft, = 2 the Wf - s are about 157 - 158 Gev and the binding energies 
~ IGev. To obtain the excitation energies uj ■:^mt = 175 Gev and the binding energies 
of the same order as we have, e.g. for strange BS, we should take rt ~ ^Jml. For 
example, for rt = 11 uj ^ 174.7 Gev and et=2 — 0.2Gei;. However, such big value of rt 
seems to be unrealistic. Anyway, it may be of interest that for very large mass of 
the "flavored" meson the scale of the binding energies of BS is connected with the 
scale of the ratio Fp/F^. 

Since the binding energies for strangeness, charm or bottom increase with 
increasing mass of the flavor for realistic values of constant F, the charmed and 
bottomed baryonic systems have more chances to be bound than strange BS. This is 
in agreement with the experimental fact that the difference of masses of A^ baryons 
and corresponding pseudoscalar mesons, Pp {k, D or b) decreases with increasing 



mass of flavor, just indicating that the binding energy of flavored quark in the 
baryon increases in comparison with that in meson with increasing mp. Nonzero 
quantum corrections to the energy of charmed (bottomed) baryonic systems are 
expected to be smaUer in comparison with strange baryonic systems, because of the 
greater mass of charmed (bottomed) quarks or mesons. 

The apparent drawback of the approach exploited in the present paper is 
that the motion of the system into the "charm" or "bottom" direction is considered 
independently from other motions. Consideration of the BS with "mixed" flavors 
is possible, in principle, but it demands more complicated treatment. 

The collective coordinates approach with the rigid or soft rotator variant of 
the model usually gives the masses of baryons considerably greater than the Kaplan- 
Klebanov-Westerberg model we used here, if the Casimir energies are not taken into 
account |Tl|, |T^ . One of the sources of this difference is the presence of the zero- 



modes contribution in the rotation energy of the order of Nc/Qp, see (11) |]T^, ^ |T^, 
which is absent in the oscillator model. As it was shown recently by Walliser for the 
B = 1 sector [|T3| this large contribution is cancelled almost completely by the 1-loop 
correction — zero-point Casimir energy which is of the same order of magnitude, 
|T^]. Anyway, since both approaches have led to similar results in the case of 



strange baryonic systems, we may expect the same for the case of charmlets and 
bottomlets, so, our results should be valid qualitatively, at least. 

The threshold for the charm production on a free nucleon is about 12Gev, 
for double charm - about 25.2Gev. However, for nuclei as a target the thresholds 
are much lower due to two-step processes with mesons in intermediate states and 
due to normal Fermi-motion of nucleons inside the target nucleus (see, e.g. [^2] ). 
Therefore, production of states with c = 1 and even c = 2 will be available on 
accelerators like future Japan Hadron Facility (energy SOGew), the subthreshold 
production of & = -1 systems on nuclei with small probability also will be possible. 
When the present paper was almost completed I have found a paper [|3| 



where the charmed few baryon systems have been considered within more conven- 
tional potential approach. The c = 1 system with 5 = 3 was found to be not bound, 
probably (the binding energy is of the order of iMev) , The 5 = 4 system was found 
to be bound with the binding energy not exceeding ^ lOMev. 

I am indebted to H. Walliser for useful discussions of the skyrmions quanti- 
zation at arbitrary Nc and Np, and for valuable remarks. I appreciate also critical 
remarks by A.Penin. 
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